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Abstract 

We adress the problem of interactions between the longitudinal velocity 
increment and the energy dissipation rate in fully developed turbulence. The 
coupling between these two quantities is experimentally investigated by the 
theory of stochastic Markovian processes. The so-called Markov analysis al- 
lows for a precise characterization of the joint statistical properties of velocity 
increment and energy dissipation. In particular, it is possible to determine 
the differential equation that governs the evolution along scales of the joint 
probability density of these two quantities. The properties of this equation 
provide interesting new insights into the coupling between energy dissipation 
and velocity incrementas leading to small scale intermittency. 

1 Introduction 

Important details of the complex statistical behaviour of fully developed turbulent 
flows are still unknown, cf. [1]. Especially the effect of small scale intermittency, 
i.e. the phenomenon of finding unexpected frequent occurences of large fluctuations 
of the local velocity on small length scales, is an open problem. 

It is commonly accepted that small scale intermittency is due to some kind of 
cascading process: Kinetic energy which is fed into the flow by external forces on 
some (large) scale L is assumed to be in an equilibrium with the energy being 
dissipated by viscosity on the smallest scale rj. Inbetween the integral length scale 
L and the dissipation scale 77, energy is continuously transported towards smaller 
scales by the decay of eddies [2, 3, 4]. 

The statistical properties of the turbulent cascade are usually characterized by 
means of the difference between the velocities at two points in space separated by 
the distance r, the so-called longitudinal velocity increment u(r): 

u(r) = e [v (x + re, t) — v (x, <)] . (1) 

The statistics of u(r) is commonly investigated by means of its moments S™{r) — 
(u(r) n ), the so-called velocity structure functions. In the framework of the cascade 
picture it is natural to assume that, for scales r within the inertial range 77 <C 
r <C L, the S™(r) are functions of the scale r and the rate of energy transfer e 
only: S™(r) — /(e, r). Assuming a constant rate of energy transfer, Kolmogorov 
derived his famous result (furtheron referred to as K41) for the structure functions: 
S% (r) oc r™/ 3 [2]. 

However, as pointed out by Landau [5], there is no reason to assume that a 
decaying eddy spreads its energy into equal parts. On the contrary, it is likely that 
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the decay of eddies is a stochastic process. As a consequence, the energy dissipation 
rate e(x, t) is a spatially distributed random variable. Accordingly experimental 
studies yield significant deviations from Kolmogorov's prediction [1]. 

The statistics of the energy dissipation rate is usually investigated by the scale- 
dependence of e r , the average of e over a ball with radius r located at x: 



Taking into account the stochastic nature of e r , Kolmogorov derived a modified 
result for the structure functions which is in better agreement with experimental 
data ([3], furtheron referred to as K62). 

The statstics of the velocity increment can also be characterized by means of 
its probability density functions (pdfs) p(u(r j). While for large scales r the pdfs 
are almost Gaussian, the phenomenon of small scale intermittency shows up in a 
stretched exponential-like shape of the pdfs on small scales expressing very high 
probabilities for large values of |u(r)|. Those deviations from the Gaussian shape are 
closely related to the deviations of the structure functions from the K41 prediction 
S£(r) oc r"/ 3 and can be attributed to the stochastic nature of the average energy 
dissipation rate, see for example e r [6, 7, 8]. 

These results clearly show that e r has an important influence on the statistics 
of the velocity increment. It would therefore be highly desirable to have an exper- 
imental tool which allows for a precise characterization of the interdependence of 
those two quantities. The aim of the present paper is to show that such a tool is 
given by the framework of stochastic Markovian processes. 

In a recent series of papers [9, 10, 11, 12, 13, 14], the Markov analysis has been 
applied separately to the velocity increment and the energy dissipation rate. The 
idea, inspired by the cascade picture, is to consider the velocity increment u(r) (or 
the energy dissipation rate e r ) as a stochastic variable which evolves in the scale r. 
If u(r) fulfills the mathematical condition for a Markov process, the evolution of u in 
r can be described by means of the Fokkcr-Planck equation, a generalized diffusion 
equation for the pdf p(u(r)) in the variables r and u. This equation is completely 
determined by two coefficients, drift and diffusion coefficient, respectively, which 
can be estimated from experimental data. The Markov analysis thus provides a 
possibility to measure the stochastic differential equations governing the evolution 
of the stochastic variable u in the scale r without incorporating any models or 
assumptions on the physics of the systems. A detailed explanation of this method 
is presented in [11]. 

The mathematics of Markov processes can be generalized to multidimensional 
stochastic variables. In the present paper, we use the mathematical theory of mul- 
tidimensional Markovian processes to establish an unified description of the joint 
statistical properties of the velocity increment and the averaged energy dissipation 
rate. Analysing experimental data, we derive a Fokker-Planck equation describ- 
ing the evolution of the joint pdf p(u(r), e r ) in the scale r. The properties of this 
equation provide interesting new insights into the interdependence of the velocity 
increment and the averaged energy dissipation rate. In particular we differentiate 
between deterministic and stochastic coupling. 

The paper is organized as follows: First, we shortly summarize the mathemat- 
ical formulation of multidimensional Markovian processes in section 2. Section 3 
briefly recalls the results of the one dimensional analysis for the velocity increment 
and the averaged energy dissipation rate, respectively. A short description of the 
experimental set-up and the data is presented in section 4 while the results of the 
two dimensional Markov analysis are given in section 5. A short summary and 
discussion of our results in section 6 will conclude the paper. 
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2 The mathematics of Markov processes 



This section gives a brief summary of the theory of multidimensional Markovian 
processes. For a detained discussion of the theorems summarized here, we refer the 
reader to standard textbooks like [15]. 

We consider the two dimensional stochastic variable q(r) which is defined as: 

1M = ( :S ) • (3) 



x(r) 

Here, x(r) is the logarithm of the energy dissipation rate e r normalised by its (scale 
independent) expectation value: x(r) =ln(e r /e). 

The stochastic process underlying the evolution of q(r) in the scale r is Marko- 
vian, if the conditional pdf p ( q(ri) | q(r 2 ), q(r 3 ), q(r]v) ) with n < r 2 < ... < 
fulfills the relation: 

P(q(n) | q0 2 ),q0 3 ), ...,q(r„) ) =p(q( r i) | q(r 2 ) ) • (4) 

P ( c l( r i) I Q( r 2), q( r 3), ) denotes the probability for finding certain values 

for u and a; at some scale ri, provided that the values of q at all larger scales 
r 2 ,T3, . . . ,tn are known. The condition (4) simply states that the transition from 
an eddy at scale r 2 characterized by u(r 2 ) and x(r 2 ) to the "state" q(ri) at scale 
n should not depend on what happens at larger scales. 

If the conditional pdfs fulfill the Markov condition (4) , any TV-point distribution 
of q can be expressed as a product of conditional pdfs: 

p(q(r-i),q(r-2),--,q(rjv)) = P(q(r-i)|q(r 2 )) x p (q(r 2 ) | q(r 3 ) ) x ... 

xp(q(nv-i)|q(r„)) xp(q(rjv))- (5) 

Equation (5) is a remarkable statement: The knowledge of the conditional pdf 
P ( q( r ) I qo^o) ) (for arbitrary scales r and r with r < r ) is sufficient to determine 
any TV-point pdf of q, i.e.: the entire information about the stochastic process is 
encoded in the conditional pdf. 

Furthermore, it is well-known that for Markov processes the evolution of the 
conditional pdf in the scale r can be described by the Kramers-Moyal expansion, a 
partial differential equation for p ( q(r) | qo(ro) ) in the variables q and r. According 
to Pawula's theorem, this expansion truncates after the second term if the fourth 
order expansion coefficient vanishes. In this case, the Kramers-Moyal expansion 
reduces to the Fokker-Planck equation 1 : 
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-r—p(q,r\q ,r ) = — ( (q,r) p(q,r \q a , 

i—1 



ro) 



t^i^i^p^r^)). (6) 



Note that we changed notation writing p{ q, r | qo , ro ) instead of p ( q(r ) | q (r ) ) . 
This notation is chosen in order to indicate that the conditional pdf is a function 
of the scale r (although r, of course, is not a stochastic variable). 

By multiplying the Fokker-Planck equation with p(qo,ro) and integrating with 
respect to qo, it can be shown that the same equation also describes the r-evolution 

^^Note that wc multiplied both sides of the Fokker-Planck equation with r (in contrast to the 
usual definition as, for example, given in [15]). The factor r on the right side of eq. (6) can be 
found in the definition (8) of the conditional moments M'*'. 
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of the pdf p(q, r). Mathematically, the drift vector and the diffusion matrix 
D^ 2 ) are defined via the limit 

^ (1) (q,r) = lim Mf } (q, r, Ar), 

A) — >0 

4 2) (q,r) = A lim o M4 2) (q,r,Ar), (7) 

where the coefficients are given by: 

M^q.r.Ar) = Jl < (gj( r - Ar) - ft (r) ) | q,r > , 

M<?>(q,r,Ar) = ^ ( (gj(r - Ar) - ft (r)) x 

x (^•(r-Ar)-( Zj -(r))|q,r). (8) 

The coefficients M^ fe ^ are nothing but conditional expectation values of the veloctiy 
increment and the energy dissipation rate, respectively, and can easily be determined 
from experimental data. One may therefore hope to find estimates for the D^-* by 
extrapolating the measured conditional moments M^ fe ' towards Ar = [11], see 
also [16, 17]. 

Alternatively, the stochastic process underlying the evolution of the variable q 
in the scale r can be described by the Langevin-equation, an ordinary stochastic 
differential equation for q(r): 

d 2 
-frH(r) = fM,r) + J29ij(d,r)rj(r) (9) 

i=i 

The components of the vector T(r) represent the stochastic influences acting on 
the process. It can be shown that a variable which is described by eq. (9) is 
Markovian, if and only if the Ti(r) are 5-correlated stochastic forces with zero 
mean. If furthermore the pdf of the stochastic forces are Gaussian, i.e. if T(r) is 
^-correlated white noise, the Kramers- Moyal expansion stops after the second term 
and the conditional pdf p{ q, r|qo, tq ) is described by the Fokker-Planck equation. 

In that case the functions f(q, r) and g(q, r) can be calculated from the drift 
vector D^ 1 ) and the diffusion matrix D^ 2 '. In Ito's formalism of stochastic calculus, 
f (q, r) and g(q, r) are given by 

/i(q,r) = ^W), 

fly (q,r) = h/V) , (10) 

\ / ij 

where VD( 2 ) is to be calculated by diagonalizing the matrix D^ 2 \ taking the square 
root of each element of the diagonalized matrix and transforming the result back 
into the original system of coordinates. 

The Langevin-equation offers an alternativ way to check the Markovian prop- 
erties of a stochastic variable. The idea which was originally proposed in [18] is to 
estimate the coefficients and D^ 2 ^ from experimental data according to equa- 
tions (8) and (7) and to calculate the functions f and g according to equation (10). 
Having determined f and g in that way, the Langevin-equation (9) can be used 
to extract T(r) from the (measured) derivatives of the qi(r). If the realizations of 
the stochastic force obtained by this method are 5-correlated with zero mean and a 
Gaussian distribution, the process is governed by the Fokker-Planck equation. We 
will use this method here instead of the one proposed in [11], since for multidimen- 
sional stochastic variables it is hardly possible to check the Markov condition (4) 
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directly by means of multiconditional pdfs. Also the numerical cost for the estima- 
tion of the coefficients D( fe ) of order three and higher grows considerably with the 
order k (growth like 2 fe ). 

3 The results of the one dimensional Markov anal- 
ysis 

Here we briefly recall the results of the one dimensional analysis for the velocity 
increment and the energy dissipation rate, respectively. Detailled discussions of the 
results for the velocity increment can be found in refs. [9, 10, 11, 19], for recent 
results on the Markov analysis of the energy dissipation rate we refer to [12, 13]. 

In the case of the one dimensional analysis of the velocity increment, the Markov 
condition (4) was found to be valid for scales r^ and differences of scales Ar = 
r. i+ i — Ti larger than a certain scale l m ar, which is of the order of magnitude of 
the Taylor microscale A. The turbulent cascade thus exhibits an elementary step 
size. This phenomenon may be seen in analogy to the mean free path of molecules 
undergoing a Brownian motion. 

Within the range of scales for which the Markovian properties are fulfilled, i.e. 
for Ar > A, the conditional moments 

M {k \u,r, Ar) = -f— ((u'(r - Ar) - u{r)f \u(r)) (11) 

of the velocity increment show a linear dependence (with small second order cor- 
rections) on Ar and can thus be extrapolated towards Ar = 0. Furthermore, it can 
be shown that the fourth order coefficient does not have an important influ- 
ence on the evolution of p(u, r) and can be neglected [11]. The pdf of the velocity 
increment is therefore governed by the Fokker-Planck equation: 

-r-^p(u,r) = -A (l>< 1 >( u ,r)p(«,r)) + ^ (p^\u,r)p{u,r)) . (12) 

Drift and diffusion coefficient turn out to be linear and quadratic functions of w, 
respectively: 

D^(u, r) = — 7(r)u, 

D (2) (u,r) = a(r) - 5{r)u + (3(r)u 2 . (13) 

When the scale r is given in units of the Taylor microscale A, the linear term 7 
of exhibits an universal dependence on scale r, independent of the Reynolds 
number [19, 20]: 

7(r) « I + 0.2v^7A. (14) 

The coefficients a(r) and S(r) are linear functions of the scale r with slopes which 
decrease with increasing Reynolds number. The quadratic term (3 shows an only 
weak dependence on r but increases significantly with Re [19, 20]. 

The analogous analysis for the energy dissipation rate [12, 13] shows that the 
pdf p(x,r) of the logarithmic energy dissipation rate x(r) is also governed by a 
Fokker-Planck equation: 

-r^-p{x,r) = - J- (£>«0r,r)f>(z,r)) + ^ (V 2) {x, r)p(x, r)) . (15) 

Again, the drift coefficient shows a linear dependence on its argument x, albeit 
with a positive slope and an additional constant term: 

D^(x,r) = F(r)+G(r)x. (16) 
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In [12, 13] the scale dependence of the coefficients F and G was expressed by means 
of the invers logarithmic scale / = In (L/r). Rewritten in terms of the linear scale 
r, the parametrizations given in [12, 13] read: 

G(r) = S + -B 1 ln(I). (17) 

The diffusion coefficient D^(x,r) can in a first order approximation taken to be 
constant in x [12] but is found to depend on the scale r [13]: 

D^(x,r)«D(r)=C ^y Cl . (18) 

With a diffusion coefficient that does not depend on x, the solutions of the Fokker- 
Planck equation (15) are Gaussian [15], i.e. the pdf of the averaged energy dissipa- 
tion rate is lognormal in agreement with Kolmogorov's [3] assumption. But it also 
follows from eqs. (15), (16), (17) and (18) that the standard deviation of this pdf is 
not described by a logarithmic dependence as assumed by Kolmogorov to provide 
scaling behaviour [13]. It can also be seen from experimental data that the constant 
value for D^ 2 \x,r) according to eq. (18) is a first order approximation; the data 
presented in [12] (fig. 2) reveal a weak but nonetheless significant dependence on x. 



4 The experimental setup 

The data set used for the analysis consists of 10 s samples of the local velocity mea- 
sured in a cryogenic axisymmetric helium gas jet at a Reynolds number of 115000. 
The measurement was done in the center of the jet at a vertical distance of 40D 
(D = 2mm is the diameter of the nozzle) from the nozzle using a specially adapted 
hotwire anemometer with a spatial resolution of 3.4^m [21, 22]. We use Taylor's hy- 
pothesis of frozen turbulence to convert time lags into spatial displacements. With 
the sampling frequency of 91.9kHz and a mean velocity of 0.62m/s, the spatial 
resolution of the measurement is 6.8/xm. Following the convention chosen in [11], 
velocities and velocity increments are given furtheron in units of <7oo . It is defined 
as Coo = \/2a, where a is the standard deviation of velocity fluctuations. For the 
data set under considerations, = 0.2m/s. 

For the integral length scale L we obtain a value of 3.5mm, the Taylor microscale 
A is 120/xm and the dissipation scale 77 was estimated to be approximately 6/um. 
For further details on the experimental setup we refer the reader to [21, 22]. 

The energy dissipation rate is estimated by its one-dimensional surrogate [1] 

For experimental reasons the derivative of the velocity field has to be approximated 
by a finite difference: 

dv( x ) _ v i x + A ) - v(x) / on x 

~x a ' (2Uj 

As the dissipative scale rj is resolved by the measurement (the resolution of the 
sensor is 3,4/im, while rj is approximately 6/im), it is in principle allowed to set 
A = 1 (in units of samples). However, from the wave number spectrum of the data 
shown in figure 1, it becomes evident that the measurement is dominated by white 
noise for small scales (below A w 5), which would lead to incorrect results for the 
energy dissipation rate if A = 1 was used to estimate dv/dx. 
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Figure 1: The wave number spectra E(k) of the original data (circles) and the 
smoothed data set (squares). The wave number k is given in units of samples. The 
straight line indicates scaling behaviour according to KJ^l: E(k) oc fc~ 5 / 3 . 



We therefore applied a digital low pass filter to the data multiplying the Fourier 
coefficients with the spectral filter function 

m = T7W (21) 

The cutoff wave number was chosen to be k n = 0.2, according to the observed 
transition to white noise for wave numbers k > 0.2 (see fig. 1 which also displays 
the power spectrum of the filtered data set). Details on the method can be found 
in [23]. 

A digital filter is of course a serious manipulation of the data and it is by no 
means obvious that the smoothed data represent the " real" velocity signal in a better 
approximation than the original data set. Therefore we compared the various results 
obtained from the original and the smoothed data. In both cases we used A = 5 
for the estimation of the energy dissipation rate. It is found that for both cases 
the coefficients D( fc ) show the same functional dependencies on their arguments u, 
x and r. Thus the coefficients calculated from the original and smoothed data set, 
respectively, are identical up to constant factor [20] . For the purpose of the analysis 
presented in this paper, these effects are of no importance. We therefore restrict 
the discussion to the results obtained from the smoothed signal. For the filtered 
data set, the mean energy dissipation rate e is 0.52m 2 /s 3 . 



5 Experimental results 

Let us start with the ^-component of the driftvector. To this end, we have to 
calculate the conditional moment 

M^\u,x,r, Ar) = (u'(r - Ar) - u(r)\u(r),x(r)) (22) 

for various values of u, x, r and Ar and try to extrapolate it towards Ar = 
according to equation (7). Figure 2 shows the coefficient at scale r = L/2 for 

u = +a QO and x = +1 as a function of Ar. Over the whole range of (differences of) 
scales < Ar < 2A, the dependence of M^p on Ar can in good approximation be 
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described by a polynomial of degree two in Ar (see [16, 17]), thus allowing for an 
extrapolation towards Ar = 0. 




Figure 2: The coefficient [u, x, r, Ar) at scale r = L/2 as a function of Ar for 
u = +(700 and x = +1 (circles). The data can be described and extrapolated by a 
polynomial of degree two in Ar (line). 

Figure 3(a) shows the result of the extrapolation for (u, x, r) at scale r = L/2 
for various values of a; as a function of the velocity increment u. The w-component of 
the drift vector shows a linear dependence on the velocity increment and only weak 
variations for different values of x. To a first order approximation, the dependence 
ofD« on x can thus be neglected and we obtain: 

D^(u,x,r) w flW(u,r) = -7(r)u. 



(23) 
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Figure 3: (a): The u-component D ( u '(u,x,r) of the drift vector at scale r = L/2 
as a function of the velocity increment u for x = — 1 (circles), x — (squares) and 
x = +1 (diamonds). 

(b): The slope -y(r) of Du (u,r) as a function of the scale r (circles). 7(r) turns 
out to be identical with the slope j(r) of the one dimensional drift coefficient (line). 

The slope j(r), which for a given scale r is obtained by averaging the results of 
the fits (23) for different x, has a value of 1.1 at the scale r = L/2. By performing 
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this procedure at several scales r we are able to specify the scale dependence of 
j(r), see fig. 3(b). ft turns out that the function 7(7-) is identical with the slope 
7(7* ) of the drift coefficient of the one dimensional Fokker-Planck equation (13) for 
p(u,r). 

The finding that 7(7-) is identical with the one dimensional coefficient 7(7-) may 
be surprising at first sight, but is a direct consequence of the fact that Du^ does 
(approximately) not depend on x. This can be seen by considering the Fokker- 
Planck equation (6) for the joint pdf p(u, x, r): 



du 2 



+ ^K>P)+| ? KV)+2^(A«P) (24) 



From the two dimensional equation (24) the equation for p{u, r) can be derived by 
integrating with respect to x: p(u,r) — J p(u,x,r)dx. Assuming that the product 
(u, x, r)p(u, x, r) vanishes for x — ► ±00 one obtains: 



d ( \ 



—r-^- [ p(u,x,r)dx 
or J 

-J^ ( D uH u , r ) J p(u,x,r)dx 
+ J~2 {^J DM{u,x,r)p{u,x,r)dx 

d 2 f 

+ du 2 ~ / D u2( u > x > r )p( x > r \ u > r ) dx P( u > r )- ( 25 ) 



£>( 2 )(u,r) 

The result of this calculation is a one dimensional Fokker-Planck equation for p(u, r) 
with a drift coefficient D^(u, r) that is identical to its two dimensional counterpart 

(2) 

The second order coefficient D„„ (u, x, r) can be estimated by an extrapolation 
of the conditional second order moments 

M£) (u, x, r, Ar) = ((«'(r - Ar) - u(r) f \u(r), x(r)) (26) 

towards Ar — > in the same way as described above for the first order coefficient 
D£\ As can be seen in figure 4, the coefficient -Dl« does not depend on the velocity 
increment u and shows an exponential dependence on x: 

D^(u,x,r) = DW(x,r) = a (r) exp ( ai (r)x) (27) 

The coefficients ao and a\ defined in eq. (27) show simple dependencies on the 
scale r (sec figure 5): a is a linear function of r, while a\ is constant: 

a (r) w 0.02(^-l), 

01 w 0.9. (28) 



Note that the linear dependence of ao on the scale r as specified in eq. (28) leads to 
(2) 

negative values for D U u on scales smaller than the Taylor microscale A. This is in 
contradiction to the definitions (7) and (8) of the diffusion matrix since, according 
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Figure 4: The uu-component Duu(u,x,r) of the diffusion matrix at scale r = L/2 as 
a function of x for u = —o^, (circles) , u = (squares) and u = +000 (diamonds). 
The coefficient does not depend on u and can be described by an exponential in x 
(the straight line represents a fit to the data for u — according to eq. (27)). 

to those definitions, the diagonal-elements of D^ 2 ^ are positive quantities. The 
finding of negative values for D^u on scales smaller than A indicates that, similar 
to the one dimensional case, the Markovian property cannot be fulfilled for such 
small scales. 




r/k r/X 



Figure 5: The coefficients do (a) and a\ (b) defined in eq. (4) as functions of the 
scale r. ao is a linear function of r, a\ is approximately constant (straight lines). 

The fact that Duu(u 7 x,r) does not depend on u is a very interesting result 
since it states that the effect of intermittency is caused by the stochastic nature of 
the energy dissipation rate. This can be seen by considering equation (25) for the 
(hypothetical) K41-case that e r is not a stochastic variable but constant. In this case 
the conditional pdf p(x, r\u, r) in equation (25) is formally given by p(x, r\u, r) — 
p{x, r) = S(x — 1) and the one dimensional coefficient can easily be calculated: 
£> (2) = / D$(x,r)p(x,r\u,r)dx = D { ul{x = l,r) = D^{r). With a diffusion 
coefficient that does not depend on u the resulting Fokker-Planck equation (25) 
for p(u,r\e r = const) can be shown to be solved by a Gaussian distribution [15]. 
The effect of intermittecy can thus be clearly traced back to the statistics of e r . 
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This result is in full agreement with the results for the conditional pdf p(u(r)\e r ) 
presented in [7, 8]. 

Before proceeding with the x-components of drift vector and diffusion matrix, 

(2) 

let us draw attention to the mixed term D ux of the diffusion matrix. This term can 
be estimated by extrapolating the conditional moments 

M<» = ^ ( («'(r - Ar) - u (r)) (x'(r Ar) - x(r)) \ u(r), x(r) ) (29) 

towards Ar — > 0. Figure 6(a) shows the coefficient Mux at scale r = L/2 for 
cxcmplarily chosen values of u and a; as a function of Ar. The coefficient exhibits 
large variations with Ar and its absolute value shows a strong decrease as Ar goes 
to zero. This may be taken as a first hint that M U 2 J vanishes in the limit Ar — > 0. 




Figure 6: (a): The conditional moment M ux (u,x,r,Ar) at scale r = L/2 as a 
function of Ar for u = —Q.Qaoo,x = —3 (squares) and u = +Q.Qa 00l x = +1 
(circles). The data can be extrapolated by polynomials of degree two (lines). 

(2) 

(b): The extrapolated coefficient D u J(u,x,r) at scale r = L/2 as a function of x 
for u = —(Too (squares), u = (circels) and u = +<7oo (diamonds). 

When the extrapolation towards Ar — > is performed by fitting polynomials of 
degree two to the data in the interval < Ar < 2 A (see fig. 6a), we obtain values 

(2) (2) 

for Dux which arc small compared to the corresponding values of M ux at finite 

(2) 

Ar. Furthermore, the coefficient D ux shows fluctuations which are of the order of 
magnitude of the values themselves (see fig. 6b) . We take this as evidence that the 
mixed coefficient D U x vanishes. 

(2) 

Given that the off-diagonal element D ux of the diffusion matrix vanishes, the 
functions g(u, x, r) in the Langevin-equation (9) can easily be calculated according 
to equation (10). In Ito's formalism, g is then simply given by: 



9uu Y ^D U ui {Jux — 9xu — 0) 9xx — \j r Dxx- (30) 

With the results obtained so far the Langevin-equation (9) takes the form: 



dr 



r 



~u{r) = 2M u + .M e xpMr))r a (r), 



~x{r) = l -D^{u,x,r) + ^ l -D x 2 x\u,x,r)Tx{r). (31) 



Even though the formulation (31) of the twodimensional Langevin-equation is yet 
incomplete, it can already be used to extract the u-component T u (r) of the stochas- 

Wr' i 



tic force from measured realizations of u{r), 4-u(r) and x(r). These calculations 
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are faciliated by the fact that the diffusion coefficient does not depend on 

the velocity increment u. It is therefore not necessary to distinguish between the 
various definitions of stochastic calculus (ltd and Stratonovich, respectively), and a 
simple Euler-scheme can be applied to calculate T u (r) (for details on the numerical 
scheme see, for example, [13]). 

The stoachstic process governing the r-evolution of the stochastic variable q(r) is 
Markovian, if the stochastic force T is 5-correlated. However, as can be seen in figure 
7, the autocorrelation function Rr(Ar) of the stochastic force T u (r) exhibits nonzero 
values up to Ar A. This clearly indicates that the Markovian properties are 
fulfilled for scales r (and differences of scales Ar) larger than the Taylor microscale 
A only. Again, we recall that this result is in agreement with the one dimensional 
analysis of the velocity increment (see [11]). 




Figure 7: The autocorrelation function Rr(Ar) of the u-component T u (r) of the 
stochastic force T(r) (open circles). Rr(Ar) exhibits finite values up to scales Ar s=y 
A. 

To complete the description of the turbulent cascade by the Markov analysis, 
we still have to determine the coefficients and D^J by extrapolating the con- 
ditional moments (u, x, r, Ar) and M^J (u, x, r, Ar) , respectively. As shown in 
fig. 8(a) for exemplarily chosen values of u, x and r, the first order coefficient M^ 1 
can be described by a polynomial of degree two in Ar over the whole range of scales 
< Ar < 2A, which again allows for an extrapolation towards Ar = 0. 

The coefficient turns out to be a linear function of the logarithmic energy 
dissipation rate x(r) and does not depend on the velocity increment u, see fig. 8(b): 

(u, x, r) = £>« (x, r) = F(r) + G(r)x (32) 

Since the coefficient does not depend on u, it has to be identical with the 

coefficient of the one dimensional Fokker-Planck equation (15) for p(x,r). 

This follows from an analogous consideration leading to equation (25). 

Plotting F and G as functions of the linear scale r, we find that their scale 
dependence is best described by (see fig. 9): 

F{r) w 0.05 - 0.04 In 

/ r n. 0.57 

G(r) w 0.03 (-) . (33) 
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Ar/X 



Figure 8: (a): The conditional moment '(u,x,r, Ar) at scaler = L/2 as a func- 
tion of Ar for u = — 0.6(Too,.t = +1 (circles) and u = +0,60-00,2; = +1 (squares). 
The data can be extrapolated by polynomials of degree two (lines), 
(b): The extrapolated coefficient D^\u,x,r) at scale r = L/2 as a function of x 
for u = —(Too (squares), u = (circles) and u = +(Too (diamonds). 

Note that these parametrizations differ from those given for the one dimensional 
coefficients in equation (17). However, the discrepancy between the results given 
in (33) and (17) must not be taken too serious. In [13] F and G were plotted 
in terms of the invcrs logarithmic length scale / = ln(L/r), which may suggest a 
different functional dependence of those coefficients on the scale than it is found 
here. Furthermore, a data set at R\ = 341 was used, whereas for the data set 
used here the Taylor-Reynolds number is 463. In addition, as mentioned above, 
numerical values such as those for F(r) and G(r) strongly depend on the method 
chosen to estimate the derivative |^ (see chapter 4). 




Figure 9: (a): The constant term F(r) of the drift coefficient D x as a function of 
the scale r. The full line indicates a fit according to eq. (33). 

(b): The slope G{r) of as a function of r . Full line: fit according to eq. (33); 
dotted line: fit according to eq. (17). 

The last coefficient which remains to be calculated is the ra-component of the 
diffusion matrix. Unfortunately, the estimation of D^J (u, x, r) from the conditional 
moment 

Mj 2 x ) («, x, r, Ar) = JL - Ar) - x(r)f \u(r), x(r)) (34) 
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yields several difficulties. When plotted as a function of Ar, the conditional moment 

( 2) 

MxJ shows a strong decrease as Ar goes to zero (see fig. 10): The extrapolated 

(2) (2) 

values of D xx are small compared to the values of M xx for finite values of Ar. 
Accordingly, the coefficient D x 2 } at scale r = L/2 docs not exhibit systematic de- 



0.05 




o.oo I — i — . — i — i — i — i — i — i — i — i — i — i — . — i — i — . — i — i — . — I 
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Figure 10: The conditional moment M xx (u, x, r, Ar) at scale r = L/2 as a function 
of Ar for u = —0. 60-00,2; = +1 (circles) and u = +Q,Qcr oc ,x = +1 (squares). The 
data were extrapolated by polynomials of degree two (lines). 

pendencies on its arguments x or u, see fig. 11(a). Although the values seem to 
decrease for large values of \x\ with a maximum at x ~ —2, the considerable scatter 
of the data also allows to assume a constant value for D X 2 J . A different result is 

obtained for smaller scales, as shown exemplarily for r = 4A in fig. 11(b). In this 

(2) 

case the coefficient D xx clearly exhibits a dependence on the energy dissipation 
rate x as well as a dependence on the velocity increment u. However, from the data 
presented in fig. 1 1 (b) it is still impossible to decide how the dependence of D X 2 X is 
to be parametrized; the data can be fitted with several functions (Gaussian as well 
as Lorentzian distributions, for example) with almost equal accuracy. The ques- 
tion of whether one of these functions is to be preferred requires further detailled 
experimental as well as theoretical investigations and is the subject of an ongoing 
study. 

6 Conclusion and Comments 

Although further investigations will be necessary to complete the two dimensional 
Markov analysis of the turbulent cascade, the results obtained so far already allow 
for several interesting statements on the joint statistical properties of the longitu- 
dinal velocity increment and the averaged energy dissipation rate. 

(2) 

A particularly remarkable result is obtained for the w-component of the 
diffusion matrix, which is found not to depend on the velocity increment. This 
means that if the averaged energy dissipation rate e r was not a stochastic variable 
but constant, the pdf of u(r) would reduce to a rather simple Gaussian distribu- 
tion. The effect of small scale intermittency can thus clearly be traced back to the 
stochastic nature of e r . 

So far, our results are in accordance with earlier experimental investigations 
[7, 8] as well as with the assumptions underlying Kolmogorov's models. Significant 
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Figure 11: The xx-component DxJ (u, x,r) of the diffusion matrix as a function 
of the logarithmic dissipation rate x. (a): At scale r = L/2 for u = (circles) 
and u = +<7oo (squares), (b): At scale r = 3A for u = (circles) and u = \o~oo 
(squares). 



new information, however, are found for the u-componcnt of the drift vector. We 
found that Du\u,x,r) does not depend on the energy dissipation rate, in a first 
order approximation. This means that this coefficient is identical with the drift 
term of the one dimensional Fokker-Planck equation for p(u,r), i.e. it reveals a 
linear dependence on u with a slope 7(7-) given by eq. (14). With respect to the 
works [9, 10, 24, 25] it remains open to explain why 7 deviates from the value 1/3. 

A further new important result is that the drift coefficient does not depend 
on the velocity increment u and is a linear function of the logarithmic energy dissi- 
pation rate x as specified by eq. (32). One thus obtains a system of two stochastic 
differential equations for the evolution of u and x in r which are linked only via 
their stochastic terms while the deterministic parts of the equations for u and x do 
not depend on the other variable. 

To summarize, the Langevin-equation for the two dimensional stochastic vari- 
able q(r) = (u(r),x(r)) reads: 

-J^-u(r) = -^7(r)u(r) + m cxp (y^W) T u (r), 



~x(r) = +-G(r)x(r) + - F(r) + J - (u, x,r) Y x {r) . (35) 
or r r V r 



According to eq. (28), the coefficient m = y in the stochastic part of the equa- 
tion for u(r) is approximately constant in r (see also fig. 5). Note also that, when 
rewritten in terms of e r , the stochastic term in the equation for u(r) exhibits a simple 
power-law dependence on the energy dissipation rate: exp(aix/2) = (e r /e) 2 ~ 
\J e r /e (see fig. 5b). 

It is interesting to compare the two dimensional equation (35) with the one 
dimensional Langevin-equations for u(r) and x(r) discussed in section 3. Those 
equations read: 



d 1 / 1 

— — u{r) = — 7(r)u(r) + \ - (a(r) — 5(r)u + [3(r)u 2 ) IYr) , 
or r V r 



~x{r) = +-G(r)x(r) + - F(r) + J - D(r) T(r) . (36) 
or r r V r 



15 



Compared to the one dimensional Langevin-equation, the stochastic part of the 
equation for u(r) is considerably simpler for the two dimensional process. Note also 
that the resulting equation for the evolution of u(r) in two dimensions is symmetric 
in u — > — u: if x was constant, the solutions of eq. (35) would be symmetric 
in u, thus leading to vanishing odd order moments. Asymmetries of the velocity 
increment, which are of importance in Kolmogorov's four-fifths law, can thus clearly 
be attributed to the influence of the energy cascade. 

While the equation for the velocity increment simplifies considerably in the two 

(2) 

dimensional formulation, the diffusion coefficient D X J (u,x,r) seems to exhibit a 
more complex dependence on its arguments (see fig. 11) than the one dimensional 

(2) 

coefficient D(r). Whether the dependence of D xx on x is in fact given by the rather 
complex form indicated in fig. 11(b) is yet an open problem which requires further 
discussions of the experiemental uncertainties. However, it is also conceivable that 
an even more complete characterization of the turbulent cascade has to include the 
transversal velocity increment as well. Such a three-dimensional analysis might 
lead to further simplifications of the diffusion terms. 

Nevertheless, the results obtained so far give reason to believe that the math- 
ematical framework of stochastic Markovian processes is a suitable tool for exper- 
imental investigations concerning the joint statistical properties of velocity incre- 
ments and the energy dissipation rate in fully developed turbulence. As it is known 
that the Fokker-Planck equation holds for the conditional pdf as well, the knowl- 
edge of this equation also provides all information about any N-scale joint pdf of x 
and u, see equations 5 and 6. Thus any general moment of u(r) a x(r')P is known. 
In this way the phenomenological Fokker-Planck equation derived from the data 
provides a closed description for all moments on all scales. 
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